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Based on first-principles calculation, it has been predicted that the magnetic anisotropy energy (MAE) in 
Co-doped ZnO (Co:ZnO) depends on electron- filling. Results show that the charge neutral Co:ZnO presents 
a "easy plane" magnetic state. While modifying the total number of electrons, the easy axis rotates from 
in-plane to out-of-plane. The alternation of the MAE is considered to be the change of the ground state of 
Co ion, resulting from the relocating of electrons on Co d-orbitals with electron-filling. 
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Magnetic anisotropy energy (MAE) is a crucial param- 
eter in high density magnetic storage^ 2 -. A large magni- 
tude of MAE will enhance the robustness of one bit to 
stabilize its magnetic spins. An efficient writing process 
requires a controllable direction of easy axis, which is de- 
termined by the sign of MAE. Co-doped ZnO (Co:ZnO), 
a possible diluted magnetic semiconductor (DMS) with 
high Curie temperature (T c ), has been demonstrated to 
be a strong single ion anisotropy material 3 . Recently, 
the dependence of magnetic anisotropy on charge carrier 
concentration controlled by an external electric field has 
been observed in DMS 4 . It offers a promising route to 
apply magnetic recording capability into information pro- 
cessing unit in current semiconductor devices. Therefore, 
it is highly desirable to understand the relation between 
MAE and carrier concentration in DMSs. 

In this work, we carry out first-principles calculation 
based on density functional theory (DFT) to obtain MAE 
in Co: ZnO. Since adding the external electric-field will 
induce the accumulation of donors (electron/hole) 5 , the 
variation of carrier concentration is treated as electron- 
filling in calculations 6 . Then, we investigate the depen- 
dence of MAE in Co:ZnO on electron-filling effect. 

A 2 x 2 x 2 supercell is created from the wurtzite primi- 
tive cell of ZnO with experimental lattice parameters and 
one Zn atom is substituted by Co (Coo.o625Zno.93750). 
With fixed lattice parameters, the internal coordinates 
are optimized by using Perdew-Burke-Ernzerhof (PBE) 7 
parameterization of generalized gradient approximation 
(GGA) as implemented in VASP package 8 . 

The electronic structure and magnetic properties are 
calculated after the structure optimization. We employ 
a corrected-band-gap scheme (CBGS) 9 , enhancing the 



TABLE I. The spin moment (Ms) and orbital moment (Ml) 
(all in /xb) and MAE (meV) with electron- filling (SN), where 
SN is the electron number added. The AMl is calculated as 
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bandgap of pristine ZnO to 2.8 eV. Besides, the Hubbard 
U (2 eV) is applied to Co ^-orbital 6 . The cut-off energy 
of plane wave is 500 eV. For the integration in Brillouin 
zone, the tetrahedron method with Blochl corrections^ 
is employed under a 5 x 5 x 3 A:-mesh. The accuracy of 
electronic iterations is up to 10 -6 eV. Then, spin-orbit 
coupling (SOC) is performed in non-self-consistent calcu- 
lations based on the wavefunctions and charge distribu- 
tion from self-consistent calculations without SOC. The 
MAE is estimated as AE = Ejioo] — £[ooi] ? where £[ioo] 
and £"[001] are the total energy in [100] and [001] magne- 
tization directions, respectively. Finally, the calculations 
above are repeated by adding extra electrons (SN) from 
-1.0 to 2.0 e with respect to the initial value. 

In Table HI we summarize the results of MAE, spin 
moment (Ms) and orbital moment (Ml) with electron- 
filling. Both of the total Ms and the Ms on Co drop 
monotonously when increasing the number of electrons 
(Fig. [Q). The MAE of charge neutral Co:ZnO is nega- 
tive, implying an "easy plane" magnetic state, which is 
in accordance with experimental results shown in Ref. S 
However, when increasing and decreasing one or more 
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FIG. 1. The dependence of (a) spin magnetic moment and (b) 
MAE on electron-filling. The positive MAE indicates that the 
easy axis parallels to the anatase c axis. 
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FIG. 2. The density of states with CBGS and U Co -d, (a) 
SN = —1, (b) SN = and (c) SN = 1. The Co d-orbitals are 
filled in red color. The vertical dash line refers to the Fermi 
level. 



electrons, the sign of MAE alters, yielding the direc- 
tion spin preferred switches from in-plane to out-of-plane. 
Meanwhile, the magnitude enlarges to more than 0.2 
meV, even up to 0.319 meV for SN = 2. Besides, it 
should be noted that the Ml in the direction of prefer- 
ential spin axis is larger than that in another, as given 
by Bruno's model 11 . To interpret the variation of MAE 
with electron-filling, the electronic structure for different 
electron number is first checked. The density of states 
are calculated and projected to Co ^-orbital for SN = 1, 
and 1. For the neutral Co:ZnO [Fig. Hfb)], the gap of 
pristine ZnO is preserved. The impurity bands, formed 
by the hybridization between the levels of the host Zn 
vacancy (oxygen dangling bonds) and Co ^-orbital, lie in 
the bandgap. The Co d-orbital in majority spin, which 
is just above the valence band maximum (VBM) of the 
ZnO host, is completely occupied. The topmost occupied 
state is Co e state in minority spin with 1.5 eV higher 
than VBM. Additionally, the e — ti split is about 2.3 
eV, and the state in minority is pushed into conduc- 
tion band (CB). The system is insulator. Our calculation 
confirms the electronic configuration of Co 2+ , <i 7 , with a 
high spin state in neutral Co: ZnO, which is in agreement 
with previous work 6 . 

When decreasing one electron (SN = — 1) [Fig. [2fa)], 
the initial double degenerated and fully occupied e state 
of Co in minority spin splits into two bands, one strongly 
hybrids with valence band of ZnO (O p-orbital) and the 
other is empty, just locating above the Fermi level. For 
SN — 1 [Fig. Efc)], doping with one electron, the empty 
ti state in minority spin is partially filled. The system 
becomes metallic. Compared to the neutral system, the 
electrons relocate in the d-orbital of Co in minority spin 
when SN = — 1 and 1, which finally results in the change 
of the ground state of Co ion. 

As proposed by van Vleck 1 — , the MAE is induced by 
the relativistic SOC interaction. Specifically, for 3d tran- 
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FIG. 3. Bandstructure of minority spin along the high sym- 
metric line in Brillouin zone, (a) SN = — 1, (b) SN — and 
(c) SN — 1. The red color of the point refers to the state 
with contribution from Co d-orbital more than 2%. The hor- 
izontal line at eV is the Fermi level, (d) Schematic diagram 
(minority spin configuration) of the perturbation determining 
MAE in Co:ZnO with electron-filling. The MAE for SN = is 
decided by the perturbation between the e level and ti level. 
However, that perturbation for SN — — 1 and 1 exists within 
e and ti manifold, respectively. 
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FIG. 4. The dependence of the total energy on the angle (0) 
between magnetization direction and the wurtzite c axis, (a) 
SN = and (b) SN = 1. 



sition metal elements, the interaction is weaker than 
the crystal-field split. Therefore, according to the sing- 
ion anisotropy theory^, MAE can be estimated by the 
second-order perturbation between the ground state and 
excited states of magnetic ion. Since the Co ion exists as 
a high spin state in Co:ZnO (Fig. [2]), which implies the 
exchange split is larger than crystal- field effect, and the 
^-orbital in majority spin is fully occupied, the variation 
of ground state of Co ion under electron-filling originates 
from the electron relocating in minority spin. 

To describe this physical scenario straightforwardly, we 
calculate the bandstructure of minority spin for different 
SN and at the same time mark the state with contribu- 
tion from Co ^-orbital more than 2%, as presented in 
Fig. [31 For SN = [Fig. [3(b)], the almost two-fold de- 
generated e level in ground state is completely occupied, 
while the level in conduction band is empty. There- 
fore, the excited states associate with electrons hopping 
from the e level and ti level, i.e., the MAE for charge neu- 
tral Co:ZnO is determined by the perturbation between e 
level and t<i level. However, as shown in Fig. [3(a) and (b), 
this perturbation for SN = — 1 and 1 exists within e and 
ti manifold, respectively. We conclude that the change 
of the perturbation process [Fig. [3(d)] finally results in 
the reorientation of easy axis with respect to SN = 0. 

Further, we calculate the total energy as a function of 
the angle (9) between the magnetization direction and 
the wurtzite c axis. Obviously, as shown in Fig. [2J we 
predict an easy-plane and easy-axis magnetic anisotropic 
state for SN = and 1, respectively. Moreover, the trend 



of the data implies that the contribution to MAE mainly 
comes from the second-order uniaxial anisotropy. By fit- 
ting the total energy to Co + Ci cos 2 curve, the standard 
error for the two coefficients, Co and Ci, up to 10 -6 sup- 
ports the former speculation. 

In conclusion, first-principles calculations are carried 
out to investigate the relationship between MAE and car- 
rier concentration in Co:ZnO. Results show that Co:ZnO 
in neutral possesses an "easy plane" magnetic state. 
While modifying the total number of electrons by one or 
more, the system alters the sign of MAE, switching into 
"easy axis" magnetic state. The electronic configurations 
on Co ion with electron-filling are checked. We conclude 
that, due to the electron relocating on Co d-orbital, the 
perturbation between the ground state and excited states 
for different total number of electrons has changed, fi- 
nally reorienting the direction of spin axis. Further, we 
investigate the dependence of total energy on the mag- 
netization direction. By fitting the data, a second-order 
uniaxial anisotropy has been found in Co:ZnO. 
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